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ABSTRACT 

We study energy dependencies of quasi-periodic oscillations from a number of black 
hole X-ray binaries. The selected sources were observed by RXTE at time periods close to 
state transitions and showed QPO in the 1-10 Hz range. We have constructed QPO root- 
mean-square energy spectra, which provide information about underlying physical process 
leading to QPO generation. These spectra show an interesting anti-correlation with the time 
averaged spectra. The QPO rm.s. spectra are harder than the time averaged spectra when the 
latter are soft, while they are softer than the time averaged spectra when the latter are hard. We 
then discuss these observational results in the context of simple spectral variability models. 
Hard QPO spectra can be produced by quasi-periodic modulations of the heating rate of the 
Comptonizing plasma, while soft QPO spectra result from modulations of the cooling rate by 
soft photons. 

Key words: accretion, accretion discs - instabilities - radiation mechanisms: thermal - bina- 
ries: close - X-rays: binaries 



■ 1 INTRODUCTION 

A common feature in X-ray power spectra from many X-ray bina- 
ries (XRB) are quasi-periodic oscillations (QPO). Many "flavours" 
of QPO have been identified, with some correlations between their 
' frequencies (Psaltis, Belloni & van der Klis 1999). For example, 
, low frequency (/ of a few Hz) QPO in black hole XRB seem to ap- 
' pear preferentially when the state of the source changes (Rutledge 
et al. 1999). They are mostly observed in luminous states, high or 
very high state, when the energy spectrum is rather complex, show- 
ing both disc thermal component and high energy Comptonized 
component (see QPO observations review in Wijnands 2001). Mod- 
els of QPO so far concentrated on identifying the frequencies of 
oscillations, with very little if any attention paid to the fact that it is 
the hard X-rays which are modulated (see QPO models review in 
Psaltis 2001). The results of Gilfanov, Revnivtsev & Molkov (2003) 
and Revnivtsev & Gilfanov (2005) are perhaps the only ones where 
the energy spectra of QPO are discussed and interpreted. These au- 
thors applied Fourier frequency resolved spectroscopy (Revnivtsev, 
Gilfanov & Churazov 1999, 2001) to show in particular that the 
kilo-Hz QPO spectra in neutron star XRB are the same as spectra 
at other Fourier frequencies. This can be interpreted as common 
spatial origin of the X-rays showing the usual broad band noise 
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variability and the more coherent quasi-periodic variability. The 
spectral shape of these components is consistent with Comptonized 
emission from a boundary layer. 

Fourier spectroscopy techniques are particularly suitable for 
analysis of QPO. However, the approach of Revnivtsev et al. (1999) 
requires integrating the power density spectra over a narrow fre- 
quency interval around a chosen frequency resulting in contamina- 
tion of the Fourier frequency resolved spectrum of the QPO by a 
contribution from the broad band variability component. Thus, in 
this paper we construct a root-inean-square spectrum, where each 
point is the r.m.s. variability integrated over the component (e.g. a 
Lorentzian) describing the QPO in the broad band power spectrum. 
Such energy dependences of QPO amplitude in XRB has been in- 
vestigated before (e.g. Belloni et al. 1997; Cui et al. 1999; Rao et al. 
2000; Tomsick & Kaaret 2001; Rodriguez et al. 2004; Casella et al. 
2004; see review in van der Klis 2006). The main conclusion from 
these studies is that the QPO amplitude increases with energy up to 
at least ~ 20 keV. This suggests immediately that the oscillations 
are connected with the high energy spectral component, rather than 
with the disk component. 

The significance of the r.m.s. energy spectrum is that it de- 
scribes the energy spectrum of the process responsible for a given 
variability component (e.g. a QPO), if the variability is generated 
by fluctuactions of the normalization of a separate spectral compo- 
nent (Gilfanov et al. 2003; Zycki & Sobolewska 2005, hereafter Pa- 
per I). If the QPO is produced by an additional modulation of one 
or more physical parameters determining the spectrum, its r.m.s. 
spectrum has less direct interpretation, but still it contains the sig- 
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Table 1. Log of RXTE observations 



Object 


MJD" 


yymmdd 


Observation ID 


Spectral state 


XTE J 1859+226 (I) 
XTE J 1859+226 (II) 
XTE Jl 859+226 (III) 
XTE Jl 859+226 (IV) 


51463.8 
51465.9 
51472.5 
51474.8 


991012 
991014 
991021 
991023 


40124-01-05-00 
40124-01-10-00 
40124-01-21-00 
40124-01-25-00 


soft state at the luminosity peak^ 


GRS 1915+105 (I) 
GRS 1915+105 (II) 
GRS 1915+105 (III) 


50422.0 
52768.6 
52731.5 


961204 
030509 
030402 


20402-01-05-00 
80127-03-01-00 
80127-01-03-00 


C state (class x)^ 
3 


XTE J 1550-564 (I) 
XTEJ1550-564 (II) 
XTE 11550-564(111) 
XTEJ1550-564(IV) 
XTEJ1550-564(V) 


51068.3 
51245.4 
51248.1 
51250.7 
51683.8 


980912 
990308 
990311 
990313 
000519 


30188-06-05-00 
40401-01-53-00 
40401-01-51-01 
40401-01-57-00 
50135-01-03-00 


soft state^ 
,, anomalous" state'' 

hard state^ 


4U 1630-47 (I) 
4U 1630-47 (II) 


50853.1 
50853.7 


980209 
980209 


30178-01-01-00 
30178-01-02-00 


hard state^ 



" Start of observation, MJD = JD - 2, 400, 000.5 
References to papers with spectral/timing analysis: 
1 Casella et al. (2004) 
2 Sobolewska & Zycki (2003) 
3 Rodriguez et al. (2004) 
■* Sobczak et al. (2000) 
5 Remillard et al. (2002), Kubota & Done (2004) 
Tomsick et al. (2001). Rodriguez et al. (2003) 
^ Dieters et al. (2000) 



natures of the modulation (Paper I). Energy dependencies of r.m.s. 
spectra were investigated by many authors with a goal of determin- 
ing how the spectral components vary with respect to each other. 
Recently, Gierlinski & Zdziarski (2005) investigated the r.m.s. en- 
ergy spectra of a number of XRBs in various spectral stated and in- 
terpreted them in terms of variations of plasma heating/cooling rate 
in a scenario where the time average energy spectrum consists of 
a disc blackbody and a hybrid (thermal/non-thermal) Comptoniza- 
tion. However, their r.m.s. spectra were averaged over all variability 
components (Fourier frequency). 

In this paper we perform a systematic study of energy de- 
pendencies of low-frequency QPO in black hole XRB in different 
spectral states. We analyse RXTE data from a number of sources 
observed at time periods close to or during state transitions, when 
the sources showed QPO in 1-10 Hz range. We selected data with 
good energy and time resolutions, enabling construction of power 
density spectra (PDS) in at least 12 energy channels covering the 
3-40 keV range. We then constructed the r.m.s. spectra of the QPO 
and compared them with the mean energy spectra. Finally, we com- 
pared the results of data analysis with the predictions of theoretical 
models. The models themselves are described in detail in Paper I. 



2 DATA SELECTION AND DESCRIPTION OF 
TIME-AVERAGED SPECTRA 

We made use of the RXTE archive data. We were particularly in- 
terested in observations during which a pronounced (rms > 10%) 
and coherent 1-10 Hz QPO were present in the power density spec- 
trum. A basic criterion was good timing and spectral resolution of 
the data. The log of all observations is presented in TableQ 

We extracted Standard 2 PCA spectra from top and mid layer 
from all available units. Standard selection criteria and dead time 



correction procedures were applied. The PCA background was es- 
timated using the pcabackest package (ver. 3.0), while the response 
matrices were made for each observation with pcarsp ver. 10.1 
tool. Systematic error of 0.5 per cent were added to each bin of 
PCA spectra, to take into account possible calibration uncertain- 
ties. HEXTE data (Archive mode. Cluster 0) were extracted using 
the same selection criteria as those for the PCA. Background was 
estimated using hxtback tool. For spectral modelling we used the 
PCA data in 3-20 keV range and HEXTE data 20-200 keV range. 

The spectral model consisted of a disc blackbody and ther- 
mal Comptonization components, modified by photoelectric ab- 
sorption. For Comptonization we use the THComp XSPEC model 
(Zdziarski, Johnson & Magdziarz 1996). This is parametrized by 
the asymptotic spectral index, F, seed photon temperature. To, and 
electron temperature, kTc- Disc blackbody (maximum) tempera- 
ture is denoted Tbb. In most cases we assume Tq — can 
be expected, if it is the disc photons which are Comptonized, but 
we will also consider cases where the temperatures are de-coupled. 
The THComp model computes also self-consistently the repro- 
cessed component consisting of the Compton reflected continuum 
and the Fe fluorescent Ka line (Zycki, Done & Smith 1999 and 
references therein). The main parameters here are the relative am- 
plitude of the component, R, ionization parameter of the reproces- 
sor, ^, and the inner disc radius used for computing the smearing of 
spectral features due to relativistic effects, Rin- 

Results of the fits to the time averaged spectra are given in 
Table |2| and described in detail in following Sections. All uncer- 
tainties in the Table are 90% confidence limits for one parameter 
of interest, i.e. Ax^ ~ 2.71. We allowed for a free normalization 
between the data from PCA and HEXTE. 
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Table 2. Fits to time averaged spectra 



diskbb/thComp 1 thComp 2 





bb V^e V ; 


p 




p 




R 


t 
s 


v2 /Hnf 


A 1 h J 1 o5y+izo (1) 
XTE J 1859+226 (11)" 


n -7n+0.09 
0.81 ±0.04 


- 




2.19 ± 0.08 
2.38 ±0.03 


441 18 
150(f) 


n 1 n+0.13 
0-19-0.10 
15+0-07 
'^•^''-0.03 


1 Qc;+4700 
180_iio 

(iAtli) X 103 


"11 C\l£iC\ 

(>\.mi 


XTE J 1859+226 (III)" 


0.82 ±0.03 






2.37 ±0.02 


150(f) 


0.16 ±0.03 


(2.5^1-6) X 103 


74.6/70 


XTE J 1859+226 (IV)" 


0.83 ±0.03 






2.34 ±0.02 


150(f) 


16+0-06 

^^•^"-0.02 


(2.9t|-0) X 103 


80.5/70 


XTE 11550-564(1)'= 




1.98 ±0.12 


7 5+1-2 
'•'^-0.9 


r, 90 + 0.10 

^■^•^-0.37 


150 (f) 


o-43ig:?5 




58.2/70 


XTE 11550-564(11)'= 


0.92 ±0.02 


2.86«-.- 


17+51 
^'-12 


1.88 ±0.20 


150(f) 


> 0.8 


(2.II*-? X 103 


63.5/69 


XTE 11550-564(111)*' 


0.66 ±0.01 


4.51 ±0.17 


100(f) 


2.19 ±0.05 


150(f) 


40+0-33 
"■^"-0.08 


(4.4li5g'5) X 103 


51.2/73 


XTEJ1550-564(IV)'' 


0.47 ± 0.03 


4.1 ±0.4 


100(f) 


2.21 ±0.04 


150(f) 


29+0-34 
"■^='-0.04 


(1.2t^-g) X 103 


74.5/74 


GRS 1915+105 (I)" 
GRS 1915+105 (II)'^ 
GRS 1915+105 (111)^= 


U.&»_o.08 

1 02+0-24 


r, 9-1 +0.07 

^-^1-0.05 
9 10+0. 16 


t; 1+0.4 
^•1-0.9 

^•1-0.6 


2.29 ±0.07 

9 9P.+0.09 
9 t:n+0.40 


40+29 
^"-20 

> 25 
150(f) 


'J-*' -0.20 

1 7+1-^^ 
^- ' -1.1 

> 0.25 


(2.0+*-?) X 10* 
(7.3+iOe) X 10* 
{3.5tl-l X 10* 


59.7/71 
69.2/62 
44.5/65 


4U 1630-47 (I)" 
4U 1630-47 (II)" 
XTEJ1550-564(V)" 


1 06+018 
1 05+0-22 

0-88_o:o9 






^•y°-0.04 

1.98 ±0.04 

1 fi77+0.019 
'-0.010 


> 45 
37lt° 

> 64 


10+0-06 
'-'•^'-'-0.02 

07+0-10 
^•^'-0.01 

0.11 ±0.03 


(4.9l4g) X 103 
(1.8t^-0) X 103 
(1.7+^-0) X 103 


80.9/73 
44.7/66 
41.8/69 



" Model: diskbb + thComp 
^ Model: thComp + thComp 
Model: diskbb + thComp + thComp 
(f) Parameter was fixed 



2.1 XTE J1859+226 

We study four observations of the X-ray nova XTE J1859+226 at 
the luminosity peak during its 1999 outburst. Detailed timing anal- 
ysis and identification of different types of QPO was presented in 
Casella et al. (2004). During the observations the source was in a 
very high state. The data are well described with a disc blackbody 
model with a temperature of ^ 0.7-0.8 keV and its thermal Comp- 
tonization on electrons with temperature fixed at kT^ = 150 keV 
(lower temperature, kTc ~ 40 keV is required for the first observa- 
tion). The Comptonized continuum has a photon index of F = 2.1- 
2.4. In all datasets a weak (amplitude R — 0.1-0.2), highly ion- 
ized (ionization parameter ^ — 1200-2700) reflection component 
is present. All but the first observation require an additional smear- 
ing, which we model as relativistic effects. 

2.2 GRS 1915+105 

We analysed three observations of GRS 1915+105. This is a very 
peculiar source showing complicated variability patterns and com- 
plex spectra. It seems to be always in a high/soft state, because of 
high accretion rate (Sobolewska & Zycki 2003; Done, Wardzinski 
& Gierlinski 2004). The spectral analysis of the three characteristic 
spectral states of this source, introduced by Belloni et al. (2000), 
was presented in Sobolewska & Zycki (2003). 

The spectrum from the first observation can be well modelled 
as a sum of disc blackbody and a thermal Comptonization compo- 
nent. In such a scenario, disc blackbody photons with a temperature 
of ~ 0.6 keV provide seed photons for Comptonization on elec- 
trons with kTc ~ 40 keV. The continuum is modified by reflection 
(7? ~ 0.5) from a highly ionized medium (^ ~ 10*), with relativis- 
tic smearing. An additional gaussian line at 6.4 keV is also needed 
to fit the residua. 

The same model does not describe well the two remaining 
datasets. Large residua in soft energy band remain and the model 
cannot fit simultaneously the high energy tail (above 50 keV). 



We try a two-Comptonization model and a model consisting of a 
disc blackbody and two Comptonizations. While the three compo- 
nent models give fits of comparable or marginally better quality 
ixl = 69.5/61 vs. 72.2/62 and 44.1/65 vs. 52.0/66, for datasets 
11 and III, respectively), the best fit values of the reflection ampli- 
tude are closer to 1 than for the two component models. We thus 
use the three component models in our further studies. Disc pho- 
tons temperature is fcTbb ~ 0.2-1 keV. Electron temperature in 
the soft Comptonization component is 4—5 keV, while in the hard 
Comptonization component it is 50-150 keV. Best fit reflection am- 
plitudes are high, R — 1.3-1.7, although with quite high uncer- 
tainties, and the reflection is highly ionized, ^ ^ 10*. Additional 
narrow gaussian line at 6.4 keV of EW« 50 eV is also required. 



2.3 XTE J1550-564 

XTE J1550-564 have been studied by many authors (e.g. Sobczak 
et al. 2000; Wilson & Done 2001). We chose four observations from 
the 1998 outburst (Cui et al. 1999; Kubota & Makishima 2004; 
Kubota & Done 2004) and one observation from the 2000 outburst 
(Tomsick, Corbel & Kaaret 2001). 

During our first observation (obs. 8 in Cui et al. 1999) the 
source already switched to a soft state. The spectrum cannot be de- 
scribed by a two-component continuum model (either a disc black- 
body and a Comptonization, or a two Comptonization model). A 
three component model provides a good fit: disc blackbody of 
Tbb ~ 0.6 keV and its Comptonization on cool (kTc ~ 8 keV) 
and hot (kTc — 150 keV) electrons. The data require a somewhat 
broadened, highly ionized (^ ~ 10"^) reflection (R ~ 0.4) and an 
additional narrow iron line at 6.4 keV. 

Our next three observations are representative to the ,, anoma- 
lous" very high state (Kubota & Makishima 2004; Kubota & Done 
2004), which the source entered after being in the standard very 
high state during the 1998/1999 outburst. A detailed timing analysis 
of these observations was presented in, e.g., Remillard et al. (2002). 
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Table 3. Analysis of r.m.s. spectra: frequency ranges for the QPO spectra and PCA configurations for high time resolution spectral data 



Dataset 


fqPO 


(Hz) 


binned mode 


event mode 


XTE J 1829+226 (I) 


2.2- 


-3.5 


B.8ms.l6A.0.35JI 


E.16ms.16B.36.1s 


XTE J 1829+226 (II) 


2.2- 


-3.5 


B.4ms.l6A.0.35JI 


E.16ms-16B.36.1s 


XTE J 1829+226 (III) 


4.3- 


-7.2 


B.8ms.l6A.0.35JI 


E.I6MS.16B .36.1s 


XTE J 1829+226 (IV) 


4.8- 


-7.2 


B.8ms.l6A.0.35JI 


E.16ms.16B.36.1s 




2.2- 


4.3 




F ^^9;/*; 37M 3(S 1 
i.j.U^Ua.JZ.iVl 'O.l a 


GRS 1915+105 (II) 


2.0- 


-3.6 


B.4ms.l6A.0.35JI 


E.500/XS.64M.36.1S 


GRS 1915+105 (III) 


1.5- 


-2.5 


B.4ms.l6A.0.35JI 


E.500/.ts.64M.36.1s 


XTEJ1550-564(I)^ 


2.0-3.1 


B.4ms.8A.0.35.H 


E.16^ts.l6B.36.1s 


XTEJ1550-564(II) 


4.9- 


-7.8 


B.4ms.8A.0.35.H 


E.16/ts.l6B.36.1s 


XTEJ1550-564(III) 


4.4-7.1 


B.4ms.8A.0.35.H 


E.16Ats.l6B.36.1s 


XTEJ1550-564(IV) 


5.5- 


-8.3 


B.4ms.8A.0.35.H 


E.16Ats.l6B.36.1s 


XTEJ1550-564 (V) 


1.9- 


-5.2 


B.8ms.l6A.0.35JI 


E.125/XS.64M.0.1S 



" Additional data point for / = 9.6-20.1 Hz is also present 
* Data for the QPO harmonic, / = 4-5.5 Hz are also used 



We find a very good description of the (III) and (IV) datasets with 
a model consisting of two Comptonizations. Additional disc black- 
body component is required only in the first anomalous state obser- 
vation (dataset II). The electron plasma temperatures are fcTc = 5- 
15 for the soft Comptonization and it was fixed at 150 keV for the 
hard Comptonization. The temperature of seed photons drops from 
feTbb = 0.9 to 0.5 keV. The data require a reflection component 
whose strength also decreases (from i? ~ 2 to 0.2). The reflecting 
medium is highly ionized (^ = 1300^000), and the data required 
also some relativistic smearing. The first and third datasets need a 
narrow Gaussian line at 6.4 keV. 

Our last observation was taken during the 2000 outburst, when 
the source was transiting to a hard state. This is observation 3 in 
Kalemci et al. (2001) (timing analysis) and Tomsick et al. (2001) 
(spectral analysis). We found that the data can be very well de- 
scribed with a weak disc blackbody (fcTbb ~ 0.9 keV) and thermal 
Comptonization (fcTo = 150 keV). The photon index of the contin- 
uum is very hard, T ~ 1.7. The data require a weak reflection com- 
ponent {R ~ 0.1) from an ionized medium (^ « 1700), with sig- 
nificant smearing, corresponding to relativistic smearing with inner 
disc radius of 8 i?g). The residua at 6.4 keV are fit with additional 
narrow gaussian line. 



2.4 4U 1630-47 

We chose two observations of 4U 1630-47 from the beginning of 
the 1998 outburst, when the source was still in a low/hard state. A 
detailed timing analysis of the outburst was performed in Dieters et 
al. (2000). In particular, a complex behaviour of QPO features dur- 
ing the transition was analysed. According to Dieters et al. (2000) 
classification, based on the power density properties, our observa- 
tions are an example of type A. We find that the energy spectra 
are well described by a model composed of a weak disc black- 
body component with a relatively high soft photons temperature, 
fcTbb ~ 1, and its Comptonization on hot electrons. The contin- 
uum slope is r ~ 1.9-2, intermediate between the hard slope of 
the last dataset of XTE J1550-564 and the soft state slopes. The 
data also require a weak reflection component (R ~ 0.1) from an 
ionized matter (^ = 2000-5000). 



3 QPO DATA ANALYSIS 

3.1 Light curves and power density spectra 

We used binned and event mode PCA data for the high time resolu- 
tion spectral analysis. PCA configurations used for each dataset are 
given in Table |3] We generated light curves in each energy chan- 
nel using standard tools from the FTOOLS package. We used the 
binned mode data in the range 3-13 keV and the event mode data 
above 13 keV. We computed the power density spectra in the Leahy 
normalization with the white noise subtracted in POWSPEC. Then, 
we corrected each PDS for background and renormalized to the 
Miyamoto normalization (i.e., (r.m.s./mean)^), by multiplying by 
^'^.P , where S and B are source and background count rates, re- 
spectively (Berger & van der Klis 1994). 



3.2 QPO r.m.s. spectra 

We modeled the PDS either as a broken power law continuum with 
a Lorentz QPO profile, or as a sum of Lorentz profiles describing 
both the continuum and the QPO (e.g. Nowak 2000). These are 
phenomenological descriptions, however, since a Lorentz function 
is a Fourier transform of a damped harmonic oscillator, it can be 
assumed (in particular in the latter case) that each Lorentzian in 
the PDS is a signature of a quasi-periodic process with a different 
degree of coherence. The QPO energy spectrum was then obtained 
by integrating the QPO Lorentzians over / for different energies. 

We constructed response matrices for binned and event mode 
data using pcarsp (the background correction was performed at 
the stage of constructing PDS) and we fit the jff-spectra in XSPEC 
in order to quantify the spectral trends. We used data in a relatively 
wide energy range, ~ 3-40 keV, but since the PCA data become 
background dominated above 20-30 keV, any residua and features 
at 20-40 keV should be treated and interpreted with caution. 



4 RESULTS 

Here we compare the energy dependent QPO r.m.s. spectra to the 
time averaged spectra. The latter were described in earlier sections 
on individual objects. The results are presented in the form of ra- 
tios of QPO data to the model components of the time averaged 
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Figure 1. A comparison of the QPO and mean energy spectra in the soft spectral state of XTE J 1859+226. Upper panels: time averaged and QPO energy 
spectra, unfolded using inodels described in Sec.|2]and Sec. 14.31 respectively. Solid curves show the total spectra, the dotted curve (magenta online) is the 
disc blackbody component, while the dot-dashed curve (blue online) is the Comptonized component. Lower panels show ratios of the QPO data to various 
model components: the total time averaged continuum (but without the reprocessed component), and to the time averaged Comptonization continuum only 
(i.e. the disc and reprocessed components removed). The QPO spectra are harder than the time averaged spectra in this soft state, and they are closer to pure 
Comptonization continua than to the total continua (i.e. the disc coinponent does not seem to contribute to QPO). The reprocessed component is present in 
two QPO spectra (I and III; see Tablel4l. 



spectra, when the model normahzation was adjusted to give small- 
est possible . These will be useful when we discuss the data in 
the context of theoretical models in the next Section. We will also 
perform formal model fits to some of the QPO spectra in order to 
more quantitatively determine their shape relative to the time aver- 
aged spectra. 

Our sample contain 14 observations of four objects. Of these 
11 observations were performed in soft spectral states while 3 in 
the hard state. 



4.1 Soft state spectra 

In the soft state of XTE J 1859-1-226 (Fig[T} the QPO spectra are 
more similar to the Comptonized component than to the total spec- 
tra. This means that the disc blackbody component is not present 
in the QPO spectra, i.e. it does not participate in the oscillations. 
The QPO spectra are harder than the Comptonized component in 
the corresponding time averaged spectrum. Ratios for dataset (III) 
clearly indicate the presence of the reprocessed component in the 
QPO spectrum, which will be discussed in more details below. 

Similarly, no disc blackbody component in the QPO spectra is 
observed for GRS 1915-F105 (Fig.|2}. In all three spectra the QPO 
data show larger deviations from the total model than from the (to- 
tal) Comptonized component. In datasets (II) and (III), where the 
time averaged spectrum requires two Comptonized components, 
the QPO are most similar to the sum of the two Comptonizations, 



rather than to either single one of them. In particular, the QPO spec- 
tra seem to be much harder than the soft Comptonization. 

We attempt to fit the QPO data with a modified time averaged 
Comptonization continuum. We first tie all the parameters of the 
QPO model, except for normalization, to the parameters of the time 
averaged spectrum. Fitting the GRS 1915-1-105 (I) dataset we also 
set the disc blackbody normalization to 0, so the only free param- 
eter is the normalization of the Comptonized component. The fit is 
very bad, xl = 267/18 dof. It improves dramatically, when seed 
photon temperature. To, is free to adjust, = 42.3/17 dof. The fit 
improves further if the continuum slope is left free, = 35.8/16. 
The resulting model has higher To « 0.95 keV compared to the 
~ 0.65 keV in the mean spectrum, and it is harder than the latter, 
r ft; 2.18 compared to ft; 2.25. The fit can be further improved if 
the reflection amplitude is free to adjust, which will be discussed 
below. 

For GRS 1915-F105 datasets (II) and (III) there is a number 
of ways the time averaged Comptonized components can be mod- 
ified to fit the QPO spectrum, since there are two Comptonized 
components in each spectrum (in addition to the disc blackbody). 
Letting To (common to the two components) to be free does pro- 
vide a good fit, with the best fit value again higher than that in 
the mean spectrum. When the seed photon temperature is fixed, 
To, QPO = To, mean, obtaining a fit of similarly good quality re- 
quires freeing at least three parameters: the two spectral slopes and 
electron temperature of the soft Comptonization. 

In the "anomalous" VHS of XTE J1550-564 the QPO spectra 
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Figure 2. A comparison of the QPO and mean energy spectra in tlie soft spectral states of GRS 1915+105 and XTE J1550-564. Panels are described in details 
in Fig.Q The dashed (green onHne) curve is the soft Comptonization component. Ratios of the QPO data to: total continuum spectrum, total (two-component) 
Comptonization, the soft Comptonization and the hard Comptonization, are plotted as labelled. 



are again harder than the time averaged spectra (Fig.|3}. The disc 
blackbody component, even if present in the mean spectrum, is not 
present in the QPO spectra. The QPO spectra slope above ~ 7 keV 
is similar to that of the hard Comptonization slope in time aver- 
aged spectrum, but there is a strong cutoff in QPO spectra below 
that energy, which formally can be described as a high value of To . 
We therefore test a model, where the QPO spectrum is equal to the 
hard Comptonization component in the mean spectrum, but their 
common seed photon temperature is untied from the disc black- 
body temperature or soft Comptonization seed photon temperature, 
To.QPO = Tchard / Tb.soft = Tbb- Such a model gives some- 
what worse fits to the combined datasets than best models to time 
averaged data alone. Pairs of values for the mean spectrum fit 
(from Section|2} and for the combined mean spectrum and the QPO 
spectrum fits, are as follows: (65.8/69 = 0.95, 90.8/77 = 1.18), 
(51.3/71 = 0.71, 65.6/81 = 0.81) and (74.5/73 = 1.02, 
85.1/82 = 1.04), for dataset (II), (III) and (IV), respectively. For 
datasets (III) and (IV) the increase of is 14.3 and 10.6 per 10 
and 9 new dof, respectively. Thus including the QPO data does not 
significantly worsens these fits. However, when the QPO fits are ex- 
amined as data/model ratios, they are approximately the same for 
all datasets. Its the larger error bars on datasets (III) and (IV) (com- 
pared to dataset I) that make the fits acceptable. In particular the 
data show pronounced low energy cutoffs below ~ 6 keV, which is 
not well modelled by the increased To if the latter is also to fit the 
mean spectrum. We thus conclude that the QPO data are unlikely 



to be described simply by the same shape as the hard Comptoniz- 
ing continuum. This implies that more complicated models, which 
include spectral variability, must be employed to explain the QPO 
in these datasets. We intend to explore such models in more details 
in a future paper. 



4.2 Hard state spectra 

In the hard spectral state of XTE J1550-564 (dataset V; Fig.|3, the 
QPO spectrum is somewhat softer than the total spectrum. When 
the model for the time averaged spectrum is fitted to the QPO 
data, the best fit model does contain the disc blackbody compo- 
nent (Tbb ~ 0.89 keV). When the disc blackbody component is 
removed and the seed photon temperature, Tq, is fixed at Tbb, the 
fit is worse by Ax^ ~ 7.2 for one more dof (x2 = 13.1/10 vs. 
5.7/9). Now, if To is free to adjust the fit improves to x2 = 5.7/9, 
but the new To = 0.54t°-j3 is lower than Tbb. Thus, contrary to 
the QPO spectra in soft states, here the QPO spectrum either con- 
tains the disc blackbody component, or the seed photon tempera- 
ture in QPO spectra is lower than that in time averaged spectra. In 
either case the QPO spectrum slope is softer than the time averaged 
spectrum. 

In the other two hard state observations, of 4U 1630-47, the 
QPO slope seems to be similar to the Comptonization slope in 
the time averaged spectrum, but the residuals suggest somewhat 
higher To (Fig.|4j. When the time averaged model is renormalized 
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Figure 3. A comparison of the QPO and mean energy spectra in the "anomalous" very high state of XTE J1550-564. Panels are described in details in Fig.lTI 
Ratios of the QPO data to: total Comptonized continuum, soft Comptonized component and hard Comptonized component are plotted as labelled. 



(r and To fixed) to the QPO data, the fits have xl = 27.5/19 
and Xv = 25.4/19 for datasets (I) and (II), respectively. Allow- 
ing r and To to adjust we obtain best fits with xt ~ 11.0/17 and 
xt = 17.5/17, and the confidence contours are plotted in Fig.|5| 
The contours suggest that the QPO spectra are again somewhat 
softer and indeed have higher To than the time averaged spectra, 
however the difference between the spectral slopes in not highly 
significant. We note then that for 4U 1630-47 the QPO spectra are 
closer to the time averged spectra than those of XTE JI550-564 
(V). Since the time averaged spectra 4U 1630-47 are softer than 
the spectra of XTE JI550-564 (V), we might be observing here a 
transition from the hard state behaviour (as in XTE J 1550-564) to 
the soft state behaviour, where the QPO spectra are harder than the 
time averaged spectra. 



4.3 The reprocessed component in the QPO spectra 

We have checked if the reprocessed component is present in the 
QPO spectra. Firstly, we assume the time averaged continuum as 
a model for the QPO, and let the reflection parameters, R, ^ and 
Rin to be free. Such models do not provide good fits, with only 
two cases when xl < 1-5 (Tabled columns a and b). Therefore, 
although the reprocessed component is present in the model, it is 
difficult to assess its statistical significance. In the two cases where 



the fits are acceptable (XTE J 1859-1-226, datasets I and IV), the 
reflection is significant. 

Secondly, we fit the QPO data with the THComp model al- 
lowing all its parameters, F, kTc and To to adjust. The continuum 
is thus described correctly, which allows for meaningful determina- 
tion of the presence and parameters of the reprocessed component. 
We then add the reprocessed component, assuming first a cold mat- 
ter (^ = 0), and then allowing ^ to adjust as well. In Table |4| we 
present the quality of the initial model (i.e. no reflection) fits, and 
the Ax^ values when the reprocessed components are added. 

Adding the reprocessed component improves the fits in most 
of the datasets, with ionized reprocessing giving better fits than cold 
reprocessing. Reduction of x^ values by 10-60 means that the re- 
processed component is highly significant. 

Reflection amplitudes in the QPO spectra are generally rather 
poorly constrained. In the three anomalous very high state spectra 
of XTE J1550-564 the best fits are obtained for a pure reflected 
component. However, the ionization parameter in the pure reflec- 
tion models is rather high, ^ ^ 10''', which means that the Fe re- 
processing features are not very prominent. The 90 per cent lower 
limits to R are rather low, nevertheless, reduction of by 11-22 
means that the component is highly significant in at least two of 
those spectra. 
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Figure 4. A comparison of the QPO and mean energy spectra in the hard spectral states. Panels are described in details in Fig.Q Ratios of the QPO data to: 
total continuum spectrum and the Comptonized continuum are plotted as labelled. 
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Figure 5. Confidence contours (Ax^ = 2.3, 4.6 and 9.6) in the T-Tq 
plane for the QPO spectra from the observations of 4U 1630-47 in the hard 
state. 'X' mark the best-fit parameters to the time averaged spectra. The 
QPO spectra are somewhat softer and have higher To than the time averaged 
spectra 

5 DISCUSSION 

We have applied techniques of Fourier spectroscopy to the analysis 
of X-ray spectra of a number of black hole X-ray binaries showing 



QPO in their PDS and we have analysed the QPO r.m.s. spectra. 
We find that in general the quasi-periodically variable part of X-ray 
spectra are consistent with thermal Comptonization and they do 
not require the disc thermal component for their description. This 
extends previous studies of Churazov et al. (2001) (for soft state 
of Cyg X-1), and it argues for the origin of the amplitude of the 
variability (including the QPO) in the hot plasma rather than in the 
cold disc. 

Comparing the time averaged and the QPO spectra we ex- 
plicitly find that only one component, a thermal Comptonization 
(with its Compton reflection) is needed to describe the latter. How- 
ever, the QPO spectra are in most cases significantly different than 
the Comptonized component in the time averaged spectra. Qualita- 
tively, we find an interesting anti-correlation: in soft spectral states 
the QPO spectra are harder than the mean spectra, while they are 
softer than the mean spectra in hard spectral states. Together with 
two intermediate cases this forms a continuous sequence. Quanti- 
tatively, the differences between the spectra can be described by a 
difference in spectral indices or seed photons temperatures, or both, 
with the statistical quality of the various descriptions being approx- 
imately the same. 

The QPO spectra do contain the reprocessed component. In 
most datasets their presence is highly statistically significant (99 
per cent or better), although the reflection amplitude is rather poorly 
constrained. Because of that, it is not possible to determine if the 
amplitude is different than that in the time averaged spectrum. 

In Paper I we introduced a phenomenological model of QPO, 
where the background radial propagation model of broad band vari- 
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Table 4. Test for the presence of the reprocessed component in the QPO spectra 



Object 


X^'Vdof 




X^'Vdof 




R" 




/?» 


XTE J 1859+226 (I) 
XTE J 1859+226 (II) 
XTE Jl 859+226 (III) 
XTE J 1859+226 (IV) 


18.7/16 
76.5/16 
24.9/14 
17.0/15 


0.47 ±0.17 
0.72 
2.2 


27.0/16 
12.4/14 
16.5/15 
9.86/15 


-3.7 
0.0 
-2.7 
-0.85 


0.3 

0+0.3 

0.50 
0.42 


-18.2 
0.0 

-11.0 
-0.9 


> 0.8^' 
0+0.4 

> 0.2 

0+3-5 
'-'■^-0.2 


XTEJ1550-564 (I) 
XTE 11550-564(11) 
XTEJ1550-564 (III) 
XTEJ1550-564 (IV) 


17.5/7 
65.1/6 
84.6/8 
18.0/6 


1.7 
0.6 
1.3 
3.0 


7.2/7 
30.8/7 
15.7/8 
10.9/8 


-5.2 
-15.3 
-4.9 
-2.5 


0.62+°-^^ 

— 0.47 
2.1 
1.3 
1.2 


-22.1 
-11.7 
-4.8 


> 0.21'^ 

> o.ie*" 

> 0.08'^ 


XTE 11550-564 (V) 


24.2/9 


0.02 


5.7/9 










4U 1630-47 (I) 
4U 1630-47 (II) 


25.7/17 
19.5/16 


0.18 

23+°-2^ 
'^•^'^-0.13 


15.4/17 
21.4/17 


-5.4 
-1.9 


77+0.73 
'J- ' ' -0.58 

0.40 


-6.1 




GRS 1915+105 (I) 
GRS 1915+105 (11) 
GRS 1915+105 (III) 


45.4/15 
31.0/16 


1.8 
1.7 


76.4/16 
27.9/16 
7.9/15 


-30. 
-7.6 


1.7 
0.97 


-65.9 
-10.8 


> 0.8 
21+°-*® 



" X^/dof from fits of the time averaged spectrum model to the QPO spectra, allowing only the parameters of the reprocessed component and overall 

normalization to be free 
* Best fit value of R from fits " 
X^^dof from fits of the thComp model to the QPO spectra (all parameters adjustable), without the reprocessed component 
Ax^ when cold (^ = 0) reprocessed component is added, with no relativistic smearing (i.e., one additional parameter) 
" Best fit value of reflection amplitude (or 95 per cent upper limit, Ax^ = 2.71, when Ax^ = in column '') 
^ Ax^ (relative to fits '^) when ionized reprocessed component is added (2 additional parameters) 
f Best fit value of R from fits ^ 
Best fit obtained for a pure reflected component {R = +oo); 95 per cent lower limit shown 
' the disc blackbody component was included in the QPO model 



ability (Kotov et al. 2001, Zycki 2003) is supplemented with quasi- 
periodic modulation of one or more of the parameters of a Comp- 
tonization spectrum. The main parameters include the plasma heat- 
ing rate, Ihcut, the cooling rate by soft disc photons, /soft (here 
I = [L / WjGT / (mc?) is the compactness parameters), and the 
amplitude of the reflected component (describing the coupling be- 
tween the heating and cooling). Our motivation for such a model is 
that irrespectively of the physical mechanism behind the QPO, the 
modulation of the emitted X-rays can only be produced by a mod- 
ulation of one or more of the parameters actually determining the 
hard X-ray spectrum. We can now discuss our observational results 
in the light of that model. 

QPO spectra harder than the mean spectra can be obtained if 
the QPO are produced by modulation of the heating rate while the 
cooling rate does not respond fully to Ihcax modulation. This pro- 
duces spectral pivoting around the low energy end of the spectrum, 
and, in consequence, r.m.s. variability increasing with energy. Such 
a model might then be applicable to the soft spectral states dis- 
cussed in this paper. The model predicts a local maximum of the 
EW of the Fe Ka line at the QPO frequency (see fig. 3 in Paper I), 
but the quality of our data is not sufficient to test this prediction. 

QPO spectra softer than the mean spectra can be produced 
by modulations driven by the cold disc. One possibility here is a 
modulation of the cooling rate. This results in the spectra pivoting 
around an energy point intermediate between the low and high en- 
ergy ends of the spectrum (see also Zdziarski et al. 2002). In the 
limited energy band, e.g. that of RXTEfPCA, it may correspond to 
the r.m.s. spectra decreasing with energy. 

The cooling rate may be modulated with or without modula- 
tion of the seed photons temperature. The difference between the 
two is that the QPO spectrum is monotonic in energy in the for- 



mer case but it has a very deep minimum (at the pivot energy) in 
the latter case (see figs. 4 and 5 in Paper I). Additionally, QPO 
harmonics appear when Thh is modulated, since the modulation of 
/soft is then not sinusoidal (even if modulation of Tbb is). However, 
the harmonics are limited to the soft disc component only, unlike in 
some of our data, for example, XTE J 1550-564 (I). 

The other possibility of producing a soft QPO spectrum is to 
modulate the amplitude of the reflected component, R, assuming 
that it also describes the feedback between illuminating hard X- 
rays and the re-emitted soft disc photons cooling the hot plasma. 
Because of the coupling, this again generates modulations of Zsoft 
and the characteristic spectral pivoting. The most characteristic fea- 
ture of this model is the strong Fe Ka line (and the entire repro- 
cessed component) at /qpo, directly resulting from modulations 
of R. Contrary to that, the model with Zsoft modulation gives either 
generally weak Fe line, or at least a minimum at /qpo- The pres- 
ence of the reprocessed component (including the Fe Ka line) in 
our QPO spectra makes modulation of R an interesting possibility 
(see also Miller & Homan 2005). 

We note that models from Paper I involving modulation of 
/soft predict very strong soft component in the QPO spectra, which 
is not observed. It is rather hard to envision a geometry where the 
modulated soft flux would enter the hot plasma, but would not be 
directed towards an observer. Modulation of R also produces some 
soft component in the QPO spectra, but it appears less prominent 
than the one from /soft modulation models. Clearly, further devel- 
opment of the models is necessary to fully describe the data. 

Another potentially useful diagnostics is the coherence func- 
tion. While variations in different energy bands are perfectly coher- 
ent when R is modulated, the coherence function shows a complex 
behaviour, with a number of minima, when /soft and/or Tbb are 
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modulated. However, distinguishing between the models based on 
such predicted observables as the EW of the Fe Kq line, the coher- 
ence function or the hard X-ray time lags would require better data 
than currently available. 

It needs to be emphasized that majority of physical models of 
QPO envision the QPO as driven by some kind of disc oscillations. 
However, our interpretation of the energy dependencies of low-/ 
QPO, at least in soft spectral states, would rather point out to mod- 
ulations of the heating rate in the hot plasma. If so, then those mod- 
els must find a way of transferring the disc oscillations to the hot 
plasma with 100 per cent efficiency, i.e. without affecting the disc 
emission. On the other hand, most of the models were formulated 
in the context of the high frequency QPO, where the energy depen- 
dencies might be different than those for the low-/ QPO. The only 
work where a physical model of the low-/ QPO was constructed is 
that by Giannios & Spruit (2004). They consider oscillations of a 
hot ion-supported accretion flow coupled to a outer cold disc. The 
coupling is realized by a number of channels: hard X-rays illumi- 
nating the cold disc, hot protons heating the cold disc, the soft disc 
photons cooling the hot flow. Spectral variability predicted by that 
model corresponds to our case of modulating the cooling rate, that 
is, the QPO spectra are softer than the time averaged spectra. This 
would then be consistent with observed QPO behaviour in the hard 
state. The overall geometry of a hot inner flow and an outer cold 
disc fits that state too. 

Considering that the usual broad band noise X-ray variability 
is driven by instabilities and flares in the hot plasma, it is obvious 
that an important piece of physical understanding is still missing. 



6 CONCLUSIONS 

• The QPO energy spectra are harder than the time average 
spectra in soft spectral states, but they are softer than the time av- 
eraged spectra in the hard state. The QPO spectra are similar in 
slope to the time averaged spectra, when the latter are intermediate 
in slope between hard and soft (F ~ 2). 

• The disc component is absent in the QPO spectra. 

• Comparison of the observational data with simple models of 
spectral variability suggests that instabilities in the hot plasma drive 
the low-/ QPO in the soft state, while cold disc oscillations drive 
the QPO in the hai'd state. 
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